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Abstract 

We numerically study the dual field theory evolution of five-dimensional asymp- 
totically anti-de Sitter solutions of supergravity that develop cosmological singu- 
larities. The dual theory is an unstable deformation of the N = 4 gauge theory 
on M x S 3 , and the big crunch singularity in the bulk occurs when a boundary 
scalar field runs to infinity. Consistent quantum evolution requires one imposes 
boundary conditions at infinity. Modeling these by a steep regularization of the 
scalar potential, we find that when an initially nearly homogeneous wavepacket 
rolls down the potential, most of the potential energy of the initial configuration 
is converted into gradient energy during the first oscillation of the field. This in- 
dicates there is no transition from a big crunch to a big bang in the bulk for dual 
boundary conditions of this kind. 

1 Introduction 

One of the main goals of quantum gravity is to advance our understanding of the nature 
of the big bang. A particularly important issue is whether the big bang represents 
the beginning of the universe, or whether semiclassical evolution essentially continues 
further back in time, possibly involving a transition from a big crunch to a big bang 

mm. 

In recent years this problem has been studied in the context of the AdS/CFT cor- 
respondence {3j, which allows one to map toy model cosmologies with anti-de Sitter 
(AdS) boundary conditions to a dual quantum field theory living on the boundary of 
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spacetime. Explicit examples of collapsing AfiS'-cosmologies were constructed in M = 8 
supergravity in four jH [5] and five [6J dimensions, using well-defined generalizations of 
the usual boundary conditions on some of the negative mass squared scalars which al- 
low smooth, asymptotically AdS initial data to evolve into a spacelike singularity that 
reaches the spacetime boundary in finite time. The properties of these solutions and 
their duals were further explored in El |9j [10] . However, a clear understanding of the 
quantum nature of cosmological singularities is still lacking. 

The dual description of the /WS'-cosmologies found in [U El E] involves field theories 
with scalar potentials which are unbounded below that drive certain boundary scalars to 
infinity in finite time. Here we concentrate on the five-dimensional solutions, where the 
dual field theory is a deformation of M = 4 Super- Yang-Mills (SYM) theory on R x S 3 
by an unbounded double trace potential —f0 2 /2, with O a trace operator quadratic in 
the adjoint Higgs scalars. This deformation is renormalizable and the coupling / that 
governs the instability is asymptotically free [5J. 

The big crunch singularity in the bulk occurs when the boundary scalar O diverges. 
Consistent quantum evolution requires that one imposes boundary conditions at infinite 
scalar field, i.e. a self-adjoint extension of the system. With a self-adjoint extension 
the boundary evolution is unitary and predicts a perfectly homogeneous wave packet 
rolls down the negative potential and bounces back. In the bulk this behavior would 
correspond to a quantum transition from a big crunch to a big bang, as envisioned e.g. 
in ekpyrotic cosmology. 

However, it has been argued [I] that the full field theory evolution is likely to be very 
different. This is because when a homogeneous wave packet rolls down, the negative 
effective mass term in the potential amplifies long wavelength quantum fluctuations 
hereby converting the potential energy of the initial configuration into gradient energy. 
In order for the boundary theory to predict a transition from a collapsing phase in the 
bulk to an approximately homogeneous, semiclassical expanding cosmology after the 
singularity there must be a parameter regime where O returns close to its original value 
after one oscillation. However, this requires the backreaction of the inhomogeneities on 
the evolution of the homogeneous wave packet to be sufficiently small. 

Here we quantify the tachyonic amplification of fluctuations by numerically evolving 
the field theory on a lattice. To study the question of a possible big crunch big bang 
transition, one is interested in the dynamics on rather short timescales and in partic- 
ular in the first oscillation of the field. Since the occupation numbers of the relevant 
perturbation modes become rapidly very large, the system behaves classically during 
these first stages of the evolution. The scalar potential is unbounded below, but the 
self-adjoint extension essentially acts as a brick wall at infinity. We therefore model this 
by adding a steep regularization term to the unbounded potential and considering the 
limit in which the UV regulator 7 introduced this way is taken to be very small. We find 
that in the 7 — > + limit, the tachyonic amplification of perturbations is a very efficient 
mechanism to produce particles and prevents the homogeneous mode from rolling back 
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up the potential. We conclude, therefore, that for 'brick-wall' boundary conditions at 
infinity and for the class of models considered here, the dual description predicts no 
transition from a big crunch to a big bang in the bulk. 

At the same time, however, our simulations provide some insight in what kind of 
models might exhibit a different behavior near the singularity. We discuss this briefly 
at the end of this paper. 



2 AdS Cosmology and its Dual Description 

2.1 AdS Cosmology 

Our starting point is J\f = 8 gauged supergravity in five dimensions, which is thought to 
be a consistent truncation of ten-dimensional type IIB supergravity on S 5 . This admits 
a consistent truncation to gravity coupled to a single S'0(5)-invariant scalar ip. The 
action then reduces to 

S = jd'x^-g Qr( 2 > - 1(^) 2 + (I5e 2 ^ + lOe"^ - e^)) (1) 

where 7 = a/2/15 and where we have chosen units in which the five-dimensional Planck 
mass is unity. The maximum of the potential at (p — corresponds to the AdS§ vacuum 
solution. Small fluctuations around this have m? = —^/R\ds^ which saturates the 
Breitenlohner-Freedman (BF) bound in five dimensions. 

In all asymptotically AdS solutions, the scalar field (p decays at large radial coordi- 
nate r as 

a(t,Q) lnr + (3(t,Q) 

where (t, Q) are the time and angular coordinates of the boundary K x S 3 of the AdS 
cylinder. To define the theory one ought to specify boundary conditions at r = 00 on the 
metric and scalar field. This amounts to specifying a relation between a and in (|2]). 
The usual boundary conditions correspond to taking a = 0, and leaving (3 unspecified. 
Alternatively, one can adopt boundary conditions of the form [TT] 

where W(/3) is an essentially arbitrary real, smooth function. Boundary conditions of the 
form (|3| are invariant under global time translations. The conserved energy associated 
with this depends on W, but is well-defined and finite. In this paper we adopt scalar 
field boundary conditions of the form a = f(3, with / an arbitrary constant. The 
corresponding asymptotic form of the metric and the expression of the conserved mass 
can be found in [121 [TBI. 



V-^^ + ^ (2) 
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For small, positive values of /, solutions were found [4J in which smooth, spherically 
symmetric initial data of approximately zero mass evolve into a big crunch - a spacelike 
singularity that reaches the boundary of AdS in finite global time. These solutions 
can be viewed as a five- dimensional version of open FLRW universes in which ip rolls 
down the negative potential, causing the scale factor a(t) to vanish in finite time. In 
terms of global time t it was found that ft ~ ft(t = 0) / (cos t/RAds) 2 - Hence ft — > oo as 
t/RAds 7 r/2, which is when the singularity hits the boundary jl]. 



2.2 Dual CFT description 

Classical Evolution 

With the usual a = boundary conditions, the dual field theory is Af = 4 super Yang- 
Mills theory. The bulk scalars that saturate the BF bound in AdS correspond in the 
gauge theory to the operators j~ Tr[$ l $ J — (l/6)5 y '<& 2 ], where are the six scalars in 
J\f = 4 super Yang-Mills. The SO (5)-invariant bulk scalar (p that we have kept in Q 
couples to the operator 

i=2 

In general, imposing nontrivial boundary conditions a{ft) in the bulk corresponds 
to adding a multi-trace interaction W(0) to the CFT action [T4"l [15], such that after 
formally replacing O by its expectation value ft one obtains (J3]). Hence the boundary 
conditions a = fft that we have adopted correspond to adding a double trace term to 
the field theory action 

S = S - W{0) = S + {- [o 2 . (5) 



O = — Tr 

N 



(4) 



The operator O has dimension two, so the extra term is marginal and preserves confor- 
mal invariance, at least classically. 

We have taken the constant / to be small and positive in the bulk. The term we 
have added to the CFT action, therefore, corresponds to a negative potential. Hence it 
is plausible that the dual field theory admits negative energy states and has a spectrum 
that is unbounded below. This would mean that the usual vacuum must be unstable, 
and that there are (non-gravitational) instantons which describe its decay. After the 
tunneling, the field rolls down the potential and becomes infinite in finite time. This 
corresponds to the big crunch singularity in the bulk and provides a qualitative dual 
explanation for the fact that the function ft of the asymptotic bulk solution (|2| diverges 
as t/R — > 7r/2, when the big crunch singularity hits the boundary. Indeed, ft is inter- 
preted as the expectation value of O in the dual CFT, and the bulk analysis predicts 
that to leading order in 1/N, (O) diverges in finite time. 
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For our purposes it suffices to concentrate on the steepest negative direction of the 
potential. Fluctuations in orthogonal directions in field space acquire a positive mass 
and are suppressed. For the 5'0(5)-invariant operator we consider, the most unstable 
direction comes from the term in (JH]). From now on we focus on the dynamics of 
$1 as it rolls along a fixed direction in SU(N), i.e. = (fi(x)U, with U a constant 

Hermitian matrix satisfying Tr U 2 = 1, so that is a canonically normalized scalar field. 
The action for this scalar is then given by 

S = [ # X y/=§ (-jW - 1r<V + W) (6) 

where A = f/N 2 , = 6/R 2 is the Ricci scalar of the 3-sphere, and the coefficient of 
the <p 2 term is fixed by the conformal symmetry. 

Although this is clearly a huge simplification of the field theory, at the classical level 
it captures the bulk behavior in a surprisingly quantitative way. In particular, it admits 
the following exact homogeneous classical solution 

cosyt/R) 

which is analogous to the background evolution of the AdS cosmology considered in 
[I]. Indeed, since 2 is identified with (3 on the bulk side, the time dependence of 
this solution agrees with the supergravity prediction^ including the fact that the field 
diverges at t/R — > ir/2. Although conformal invariance is broken quantum mechanically, 
it has been shown that quantum corrections do not turn around the potential [6]. 



Quantum Evolution 

We have seen that evolution ends in finite time in the classical description of a simplified 
version of the dual field theory, in agreement with the supergravity result. It has been 
noted, however, that this conclusion changes dramatically if one considers the quantum 
mechanics of the dual model ([6]). That is, one again concentrates on the homogeneous 
mode (p(t) = x(t) only, but one now treats this quantum mechanically. 

In quantum mechanics, to ensure unitarity one constructs a self-adjoint extension 
of the Hamiltonian. This is done by carefully specifying its domain [TBI [T7] . One finds 
the center of a wave packet follows essentially the classical trajectory and still reaches 
infinity in finite time, but it bounces off infinity and reappears as a left-moving wave 
packet [IT]. Hence a quantum mechanical treatment of the homogeneous mode of the 
dual model appears to suggest that there is an immediate transition from a big crunch 
to a big bang in the bulk. 

1 lt follows from Fig 14 (b) in [6] that for sufficiently small values of / the supergravity analysis 
predicts /3 ~ 1/A. The field theory solution (l7|) reproduces this scaling. 
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However, it has been argued jl] that the full field theory evolution is likely to be 
very different. Indeed, a quantum mechanical analysis obviously neglects the possibility 
of particle production. Calculations in similar scalar field theories indicate that this 
significantly affects the evolution of the homogeneous mode. In particular, in many 
theories where the scalar field rolls down from the top of its effective potential, particles 
are produced in great numbers while the field is rolling down. This phenomenon is often 
called tachyonic preheating [181 HH]- It happens essentially because the negative effective 
mass term in the potential amplifies long wavelength quantum fluctuations. 

Tachyonic preheating is so efficient that in some theories most of the initial potential 
energy density is converted into gradient energy well before the field reaches the true 
minimum. In theories of this kind a prolonged stage of oscillations of the homogeneous 
component of the scalar field around the true minimum of the potential does not exist 
in spontaneous symmetry breaking. 

One expects that the tachyonic amplification of fluctuations similarly affects the 
evolution of the dual field theory description of the AdS cosmologies, in which the 
supergravity initial data correspond to a homogeneous field theory configuration high 
up the potential. In this case the scalar potential does not have a global minimum. 
However, the quantum mechanical self-adjoint extension - assuming for now it can 
be implemented in field theory - essentially acts as a brick wall at infinity. From a 
cosmological perspective, one is interested in the dynamics at intermediate times, namely 
whether a wave packet rolling down the potential bounces back one or several times 
before the system settles down. To answer this one must quantify the effect of the 
tachyonic instability, and the non- linear growth of fluctuations that comes with it, on 
the background evolution. This is the subject of the remainder of this paper. 



3 Dual Model 

3.1 'Brick Wall' Regularization 

We have seen that a self-adjoint extension in a quantum mechanical treatment of ^ 
acts as a brick wall at infinity. Perhaps the simplest way to try to implement this at 
the field theory level is to add a steep regularization term to the potentiaQ 

V(4>) = \<t? - \\<P 4 + e0 6 (8) 

and to consider the dynamics in the limit e — > + . It is useful to introduce the parameter 
7 = 32e/A 2 , which can be thought of as the square of a UV regularization length. In 

2 From here onwards we set the AdS radius equal to one, so that the dual field theory lives on S 3 x R 
where the sphere also has unit radius. 
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terms of 7 the potential ([8]) can be written as 

where Vi is defined as V with A = 1. For 7 < 1, the potential has two positive roots (f>± 
(see Fig[T|. The global minimum of V on the positive axis is located at <p m , and we will 
write V m = V(<p m ). 

For 7 + , to which we will refer as the scaling limit, <p_ is insensitive to the 
regularization and tends to a/2/A. On the contrary, the large (j) parameters m , (j) + and 
V m follow the scaling given by 7: 

A#~*y-\ A^^ 7 -\ W m ^-f 7 r 2 . (9) 




Figure 1: Regularized scalar potential in the 
boundary model with A = 1 and 7 = 0.4. 

In order to model the background evolution of the AdS cosmologies we consider 
homogeneous field theory solutions (pit) that start at rest at = </>_. If one neglects 
fluctuations, the homogeneous field oscillates back and forward between the positive 
roots of the potential, just as in the quantum mechanical treatment of the unbounded 
potential discussed above. The period T of the homogenous field oscillations is inde- 
pendent of the potential parameters and given by 

Hence in the limit 7 — > + the field reaches infinity exactly at t = tt/2, which corresponds 
to the time at which the big crunch singularity in the bulk reaches the boundary. This 
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suggests that for small 7 the </> 6 regularization acts very much as a brick wall. Indeed, 
the interval of time T + during which > <ft m in one oscillation scales as y/j, 

■/*,„ V^vW) V2 1 > 

However, the regularization of the boundary potential affects the evolution in the 
bulk, because it changes the bulk boundary conditions on the timelike boundary Ac- 
cording to the general prescription ^ these now read 

a = A/3 - 3e/3 2 . (12) 

For small e this can only significantly affect the evolution near where the singularity 
hits the boundary, where is large. Nevertheless this suffices to turn the cosmological 
singularity into a large, stable black hole with scalar hair [3], where the bulk scalar field 
turned on is dual to the operator Q in the boundary theory. These hairy black holes, 
which do not exist for a = A/3 boundary conditions, have been interpreted in the dual 
theory as thermal excitations about the global negative minimum V m [5j. Hence they 
are the natural end state of evolution in the bulk corresponding to wave packets rolling 
down a regularized potential ([8]). In the limit e — > where the global minimum goes to 
minus infinity the hairy black holes become infinitely large and one recovers the original 
cosmological solutions [5]. 

Hence in the black hole context there is convincing evidence that the dual system 
will eventually thermalize. However here we are more interested in the dynamics at 
intermediate times. Indeed, whether the dual theory predicts a transition from a big 
crunch to a big bang depends on whether the backreaction of the inhomogeneities pre- 
vents an initially homogeneous wave packet rolling down the potential from returning 
close to its original value after the first oscillation. 

3.2 Initial Conditions of Perturbations 



Following [TS1 [20] we set the fluctuations initially in their instantaneous adiabatic vac- 
uum, specified by 

(ty(£)ty(#)> = ^-/{k + k') 

ZU!o{K) 

5<p(k) = ±ioo (k)5(j)(k) (13) 

Because of the unstable — </> 4 term in the potential, one has tUq = V^,(0_) < 0. We 
therefore set 

^{kf={ , 2 f 2 ^' m °! (14) 

y ' 1 At + rrig k > \m \ v ' 
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to avoid that uo(k) is imaginary for k < \mo\- The initialization of the frequency (14) 
for the modes with k < \m \ is to some extent arbitrary. However, in the limit 7—7-0 
these modes, which are tachyonic at t — 0, constitute only a small part of the ensemble 
of modes that become tachyonic at some point during the evolution. As part of our 
numerical study, we have verified that changing the initialization of the k < \mo\ modes 
does not significantly affect the evolution. 

We are initializing random classical fluctuations to reproduce the quantum vacuum 
fluctuations. The consistency of this approach relies on the fact that, as we will show, in 
the subsequent evolution the occupation numbers become rapidly very large so that the 
quantum evolution behaves classically, at least during the first stages of the evolution. 

Finally we note that for fixed 7, changing A is equivalent to rescaling the width of 
the initial wavepacket. Indeed, defining X = the field equation for X is given by 

X-A X +^=0 (15) 
dx 

which is the field equation for with A = 1. The homogeneous mode also takes the 
correct initial value, since Vi(vA0_) = 0. By contrast the amplitude of the fluctuations 
5x is rescaled by a factor vA with respect to the A = 1 model (with H — c — 1): 

(S X (k) 5 X (k')) = \^~J\k + k'). (16) 

Hence to study the evolution for a range of values of A one can equivalently change the 
width of the initial wavepacket, and then go back to the physical variable 0. 

3.3 Tachyonic instability 

The effective mass m 2 = V"(<p) of is negative in a range of field values of width 
A0 ~ l/y/Xj, extending roughly from = 0_ ~ a/2/A to ~ y/3/}rf. The lowest 
value of m 2 is 

m - = 1 -^ * (17) 
This value is A-independent, as can be simply deduced from the fact that = Vi :XX . 
Modes with momenta k 2 between and k%, = —vn} min have an imaginary frequency at 
some point during the evolution and hence are tachyonic. Decreasing 7 not only widens 
the tachyonic band but also amplifies more strongly the modes within the tachyonic 
range x = k/kx < 1. This provides a mechanism for converting potential energy into 
particles, acting as a friction term for the motion of the homogeneous component of the 
scalar field. 



Eq.(17) shows that 7 acts as a UV regularization length squared. Indeed, the tachy- 
onic instability produces particles with momenta 

k < k T ~ r 1 ' 2 
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irrespective of the value of A. However, the details of the spectrum will generally be A 
dependent at times sufficiently large for the non-linearities to play a significant role. 

Systems of this kind have been examined in the context of cosmological preheating 
and dynamical symmetry breaking [181 EQj EE]. The general conclusion is that the 
tachyonic amplification is a very efficient mechanism to convert potential energy into 
particles so that it typically takes only a few oscillations before the homogeneous mode 
stabilizes around the global minimum. A sample evolution of the homogeneous field 
in our regularized model ^ of the dual theory is given in Fig [2j for A/4 = 1CT 6 and 
7 = 0.3. One sees that with these parameter values the field returns twice to close 
its original value before backreaction suddenly kicks in and dramatically changes the 
evolution. Once the homogeneous mode, which is obtained by averaging the field value 
over the lattice sites, has settled around the global minimum the system thermalizes on 
a much longer timescale through the gradual transfer of energy into higher harmonic 
modes. The central question we wish to address here is whether the first oscillation 
survives in the 7 — > + limit. 

To answer this one needs to take in account the evolution of the fluctuations at the 
non-linear level. This can be seen as follows. The linearized perturbation equation reads 

<P k + (k 2 + m 2 (t))(P k = (18) 

where m 2 (t) = V^,(</>(t)). In the absence of the regularization term, one finds that 
m 2 (t) — —6/r 2 where r = 7r/2 — t is the time away from the singularity. Hence for 



t <C 1/k, Eq.(18) becomes 



T 



(19) 



which is solved by 



= (20, 



The IR divergence in (20) obviously comes from using the near-singularity behaviour of 
the classical solution and has no physical meaning. It can be dealt with by imposing an 
IR cutoff km which mimics the actual solution. At a given r, one has 

i 1/ \,,2(\Mt = o)\ 2 ) _ , x2I4 



(5tf)(?) > 4vr / dkk "™ " 1 ~ r- 4 ~ A> 4 (21) 
J km \ KT ) 

where we have restricted the range of integration to the modes that are tachyonic at 
time r. This indicates that while the field rolls down, non-linearities become important 
when (f) ~ 1/A. 
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Figure 2: Evolution of the homogeneous component calculated by averaging 
the field value over the lattice sites, for A/4 = 10~ 6 and 7 = 0.3, taking in account 
the growth of fluctuations 5\- The colored band represents the x =t \J ($X 2 ) zone. 

4 Lattice Evolution 



To study the dynamics in the presence of the brick wall regularization, we have de- 



veloped a numerical code that solves the non-linear Klein-Gordon equation (15) on a 



3-dimensional cubical lattice with periodic boundary conditions. We employ the same 
second order algorithm used in (2TJ. Significant improvements in the numerical perfor- 
mance were made with respect to previous codes. In this section we first describe the 
numerical setup in more detail and then discuss our results for the classical evolution 
on the lattice. 



4.1 Numerical Setup 
IR cutoff 

The scalar theory ^ is denned on IR x S 3 . The conformal coupling to the curvature 
of the 3-sphere is crucial to determine the initial conditions for the evolution. On the 
other hand, in the scaling regime, the tachyonic scale l\ = k^ 2 oc 7 becomes much 
smaller than R. Hence one expects the IR cutoff given by the radius of the sphere to 
be unimportant in the 7 — )• limit, so that one can safely work on a toroidal lattice. 
To verify this is indeed the case we have compared the evolution with the lattice size 
adjusted so that the momentum gap between the homogeneous mode and the lowest 
/c-mode is the same as on the sphere, to the evolution on a very large lattice. It was 
found that in the scaling regime, the finite volume of the 3-sphere has no effectj^j 

3 We have also verified that the energy density correlation length (see [55]) is much smaller than the 
IR cutoff. 



11 



Finally we have verified that, again in the scaling regime, the statistical dependence 
on the (pseudo-) random initial conditions is negligible if one considers averaged quan- 
tities. This provides a consistency check for the semi-classical approach we employ. 

UV convergence 

We are particularly interested in the early stages of the evolution, such as the dependence 
of the first oscillation of the field on the parameter values. It turns out that, as expected, 
in this regime only tachyonic modes with x < 1 are important. This is illustrated in Fig [4] 
(right panel) where we plot the mode occupation number as a function of x after one 
oscillation. One sees the occupation numbers fall sharply for x > 1. Hence one expects 
that a lattice cutoff at xjjv — 2 should be sufficient to ensure the UV convergence of the 
lattice evolution. This turns out to be the case: in Fig [3] we show that the homogeneous 
field displacement after one oscillation becomes cutoff independent for x > 2. 

When the first oscillation ends, the spectrum is essentially cut at x ~ 1. However at 
later times, when the field relaxes around the minimum, the gradient energy is trans- 
ferred to higher x modes. This is also shown in Fig|4| right panel. This is a slow process 
of which only the first stages are given by classical evolution. However, whereas this 
is important to understand the ultimate equilibrium state of the field theory, it is not 
relevant for the 'cosmological' question whether there is a bounce in the 7 — > limit. 
For this, the field displacement Ax, shown in Fig [5] (left panel), turns out to be a useful 
variable. 

A* 

2.5 [ 

2.0 * 
• 

1.5 • 
1.0 • 
0.5 * 

0.0 0.5 1.0 1.5 2 .0^ utoff 

Figure 3: The 'homogeneous' field displacement after the first oscillation, as a 
function of the lattice UV cutoff, for A = 4 • 10~ 3 , 7 = 1(T 3 . 

4.2 Numerical results 

To evaluate the particle production and in particular its backreaction, we adopt as an 
indicator the field displacement A</> = ft 1 ' — 0_ between the homogeneous field value 
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Figure 4: The homogeneous component x as a function of time (left), and the mode occu- 
pation numbers as a function of x = k/kx, at t ~ 1.6 when the first oscillation is completed 
(dashed) and at t = 2 (solid). Here A/4 = 10" 4 , 7 = 1(T 4 . 

after the first oscillation, <jfi~\ and its initial value 0_ (see Fig[4]). By homogenous field we 
mean the average of the field value over the lattice sites. The growth of inhomogeneities 
can be quantified by evaluating the energy density p g injected in spatial gradients after 
a single oscillation. These two quantities should naturally be compared to (0 m — </>_) 
and \V m \, which determine their upper bounds in the presence of the regularization. 

To compare theories with different values of A it is useful to work with the rescaled 
quantities A% = \/X A0 and p x = Xp g . The advantage of these variables is that the 
reference quantities (% m — X-) an d Kn,i are A-independent. 

Our study covers a wide range of parameters, namely from A ~ 1CT 1 to A ~ 10~ 10 , 
and from 7 ~ 10 _1 to 7 ~ 1CT 6 . The lower bounds are imposed by numerical limitations: 
small values of A require a higher numerical precision to resolve the fluctuations, whereas 
smaller values of 7 lead to a more violent and "tachyonic" first oscillation that requires 
a better discretization. 

AdS crunch limit: 7 dependence 

For large values of 7 near 1, simulations have shown that a few oscillations occur before 
the homogeneous field reaches its equilibrium value around the true minimum of the 
potential (see Fig [2]). At first the fluctuations are in the linear regime, in which their 
amplitude is essentially multiplied by a constant (which depends on the wavenumber) 
at each oscillation. However, since the growth of perturbations is exponential, at some 
point the backreaction becomes important, rapidly dampening the oscillation amplitude. 

When one decreases 7, the potential becomes steeper and deeper. This means that 
more modes become tachyonic during the first stages of the evolution, and further that 
the low k modes grow faster. We find the energy converted in particles is not restored 
back in the homogeneous modes when the field rolls up again. As a consequence the 
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Figure 5: Homogeneous field displacement (left) and amount of gradient energy density 
(right) after one oscillation, as functions of 7 _1 . The black solid lines represent respectively 
{Xm ~ X-) and V m ,v From top to bottom, A/4 = KT 1 , 1(T 2 , 1(T 3 , 10~ 4 , 10" 6 . 



displacement of the homogeneous field after one oscillation increases when 7 decreases. 

Since for sufficiently small values of 7 the only relevant scale in the problem is set 
by 4> m , one expects that A0 oc m and p g oc V m . This is confirmed by the numerical 
simulations, as shown in Fig [5J Hence A0 as well as p g follow essentially the scaling 
that is imposed by the brick wall. The scaling limit is manifest for A0 > 0_. In this 
regime one has 

A0cx 7 -5, p 9 cx 7 - 2 (22) 



The proportionality constants in (22) depend on A and will be discussed below. 

The subsequent evolution depends on the details of the model. For intermediate 
values of 7 it takes some time for to stabilize around <p m . This process involves a slow 
cascading effect whereby the energy is transferred from the tachyonic band to a wider 
band of modes extending up to wavenumbers of a few times kx- 

However for small values of 7 the evolution proceeds in a much more turbulent 
manner. In this regime the non-linearities are important and even at times t < T/2 
we find patches with ~ m form. The inhomogeneities behave in a local manner 
till t ~ T/2, in line with the findings in [6], but for t > T/2 the evolution proceeds 
via the expansion and collision of approximately homogeneous bubbles (see also [8]). 
This produces a violent backreaction on <ft which transfers energy rapidly to high k 
modes. Consequently, no further large oscillations in occur. The bubble dynamics 
for t > T/2 is heavily dependent on the regularization, and therefore on the kind of 
boundary conditions we impose at infinity. 

The simulations have also shown that, in the small 7 regime, a certain amount of 
tunneling takes place to the true minimum of the potential located at negative <fi. In fact, 
for small 7 the potential barrier is small compared to the scale of the true minimum 
of the potential, so the potential is essentially flat around = 0. It is therefore not 
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Figure 6: Homogeneous field displacement (left) and amount of gradient energy den- 
sity injected (right) after one oscillation, as functions of A. From top to bottom, 7 = 
1(T 6 , 1(T 5 , 1(T 4 , 1(T 3 , 1(T 2 , 5 • 1(T 2 . The curves with the two largest values of 7 exhibit 
the linear regime for small A. The straight dashed lines are, respectively, the \/A and constant 
fits to the linear analysis discussed in the text. For smaller values of 7 the backreaction is 



always large and one obtains the scaling given in Eq ( 25 ) . 



surprising that local energy overdensities let bubbles where the field is negative form. 
We have checked carefully that this phenomenon is not due to numerical violations of 
energy conservation. As a matter of fact, this feature of the evolution does not seem 
crucial for our study. Indeed, an infinite potential barrier at <fi = does not significantly 



modify the scaling behavior of (22) 



Non-linearities: A dependence 

For A 7 and fixed 7 < 1, the fluctuations are small compared to <f) and behave linearly 
at least during the first oscillation^} At the linear order, both the initial conditions and 
the evolution of the <fi inhomogeneities are A - independent. Consequently we expect no 
A-dependence for p g in this regime. One should also have 



A0 



Pg 



V, 



oc v A 



(23) 



By contrast, for larger values of A the first oscillation may already involve significant 
backreaction, with A0 ^ </>_. The homogeneous field displacement is then related to p g 

as 

P 9 ~ A(A0) 4 . (24) 



4 We emphasize that, since the tachyonic instability is an extremely efficient mechanism to produce 
particles, this regime typically corresponds to values of A that are very small compared to (an appro- 
priate power of) 7. For this reason the regime in which the perturbations behave linearly for the entire 
duration of the first oscillation has turned out to be numerically inaccessible for 7 < 10~ 3 . 
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Numerical simulations have confirmed this qualitative behavior. In Fig [6] we show that 
for A <C 7 the gradient energy p g becomes independent of A and A(f) scales as \/A, as 



expected from the linear analysis (23). For larger ratios of A/7 both the homogeneous 
field displacement and the energy density in fluctuations follow a rather precise scaling 
with A which appears to be universal for sufficiently small 7. The equal spacing between 



the data sets in Fig |6| corresponds, of course, to the scaling with 7 given in (22). 

To summarize, in the scaling limit 7 — > + and well beyond the linear regime, the 
quantities of interest behave as 

X a X 13 
A0 oc — , p g oc — (25) 

72 7 



with a ~ (3 ~ —0.33, which agrees with (24) to a good accuracy. These exponents 



are non-trivial and the result of a highly non-linear evolution involving an interplay 
between tachyonic growth of individual modes and interactions. 

The 'cosmological' limit corresponds to taking 7 — > + for fixed A. In this limit, it 



follows from (25) that A(p^ and p g both diverge. Consequently, there is no meaningful 
first oscillation back to a finite value of the field, and therefore no evidence in our dual 
model for a transition in the bulk from a big crunch to a big bang. 

Potential renormalization 

A one-loop computation shows that the coupling A of the unstable — 4 potential runs 
logarithmically. For large field values the Coleman- Weinberg potential is given by 

Vr ~ -J7^<P 4 (26) 

where M is the renormalization scale. The logarithmic running effectively decreases the 
parametric divergence of the tachyonic scale with 7. At leading order in the logarithmic 
corrections one finds i 

4 oc (27) 
log 7 

Modeling the self-adjoint extension again by the addition of a 6 regularization term we 
obtain, as a toy model, 



^ ) = ^- ,,„„.^/ + ^ a (28) 



l,o A 1 

21og^ M2 

where 0o is added to produce an essentially constant A for small values of 0. 

Intuitively, one expects that the running of A leads to logarithmic corrections to 



the scalings given in (25). A numerical study of (28) confirms this: the effect of the 



16 



logarithmic running appears to add only finite 7 corrections to the scalings. This means 
the energy in particles is still determined by the tachyonic scale, and therefore our 
conclusions remain unchanged. 

5 Conclusion 

Using the AdS/CFT correspondence, the approach to a big crunch singularity in five- 
dimensional asymptotically anti-de Sitter spacetimes can be described in a dual field 
theory by an initially nearly homogeneous wave packet that rolls down an unbounded 
potential direction of a double trace deformation —f'0 2 /2 of M = 4 Super- Yang-Mills 
(SYM) theory on IxS 3 , with O a trace operator quadratic in the adjoint Higgs scalars. 
The big crunch singularity in the bulk occurs when the boundary scalar runs to infinity 
in finite time. 

Concentrating on the steepest negative direction of the potential, we have numer- 
ically solved for the dual evolution. Consistent quantum evolution requires that one 
considers a self-adjoint extension of the system. Since this essentially acts as a brick 
wall at infinity, we have modeled this by adding a steep regularization term to the un- 
bounded potential. This introduces a UV regulator 7, and the cosmological limit in the 
bulk corresponds to taking 7 — > + . 

While the wave packet rolls down the potential, the negative effective mass of the 
potential amplifies quantum perturbations which therefore behave classically, at least 
during the first stages of the evolution. We find that for small 7 this converts most 
of the potential energy of the initial configuration into gradient energy during the first 
oscillation of the wave packet. In particular, for small 7 the system tends to a universal 
'scaling limit' in which the mean field displacement A0 and the gradient energy p g after 
one oscillation scale as (an appropriate power of) the UV regulator itself. This means 
that, in the cosmological 7 — > + limit, the boundary scalar does not return close to its 
original value after the first oscillation. Instead the gradient energy diverges, indicating 
there is no transition in the bulk from a big crunch to a big bang in the models considered 
here. In the boundary theory, our simulations show a sudden phase of bubble formation 
at large field, when some regions of space bounce back while others are still going down. 

A different way to implement self-adjoint boundary conditions in field theory, based 
on an extension of the domain of the scalars involved, has been explored in [S] and 
more recently in [23]. There is some evidence [24] this allows for a smoother dynamical 
evolution across the singularity and in particular that this avoids the breakdown of 
the local behavior of the perturbations near the singularity. Hence a violent phase of 
bubble dynamics like what we observed in our simulations may not occur with boundary 
conditions of this kind. We note, however, that in the models we have considered we find 
that, in the scaling regime, perturbations are in the non-linear regime already while the 
wave packet rolls down. It seems plausible that this is largely independent of the nature 
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of the boundary conditions at infinity. One will therefore have to take in account non- 
linearities in order to determine conclusively whether the boundary scalar can return 
close to its original value after the first oscillation with different boundary conditions at 
infinity. 

A possible exception are models where the coupling A governing the potential insta- 
bility is significantly smaller than the UV regularization scale 7 for large field values. 
In this case our simulations indicate that the perturbations remain in the linear regime 
during the entire first oscillation. The backreaction of particle production can then be 
limited, at least in a certain range of parameter values. In the models discussed here, 
we have seen it is not possible to maintain A < 7 in the cosmological regime involving 
small 7. Indeed, although the coupling of the unstable double trace deformation runs 
logarithmically and is asymptotically free, we find that in order to limit the backreaction 
A would have to go to zero in the UV faster than this. Whether this can be realized in 
a quantum treatment of the three-dimensional models discussed in [22] remains to be 
seen. 

Even though in the regularized models, the expectation value of O rapidly stabilizes 
around the global minimum of the potential, the spectrum of produced particles at early 
times is essentially limited to the low tachyonic band of momenta. Subsequently the 
system thermalizes, due to the self-interaction of the field. However this takes place on 
much longer time scales and cannot be described by the classical simulations discussed 
here. The approach to a thermal final state in the boundary theory nevertheless pro- 
vides an interesting testing ground of the AdS/CFT correspondence in the presence of 
multitrace deformations. Some aspects of this will be discussed elsewhere |26j . 
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